Just how different are SU{2) and SU{3) Landau-gauge propagators in the IR regime? 
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The infrared behavior of gluon and ghost propagators in Yang-Mills theories is of central impor- 
tance for understanding quark and gluon confinement in QCD. While simulations of pure SU{3) 
gauge theory correspond to the physical case in the limit of infinite quark mass, the 5(7(2) case (i.e. 
pure two-color QCD) is usually employed as a simplification, in the hope that qualitative features 
be the same as for the SU{3) case. Here we carry out the first comparative study of lattice (Landau) 
propagators for these two gauge groups. Our data were especially produced with equivalent lattice 
parameters in order to allow a careful comparison of the two cases. We find very good agreement 
between SU{2) and 5'!7(3) propagators, showing that in the IR limit the equivalence of the two 
cases is quantitative, at least down to about 1 GeV. Our results suggest that the infrared behavior 
of these propagators is independent of the gauge group SU{Nc), as predicted by Schwinger-Dyson 
equations. 



PACS numbers: 11.15.Ha 12.38.Aw 14.70Dj 



INTRODUCTION AND MOTIVATION 



Despite recent progress, the infrared structure of Yang- 
Mills theory is still not fully understood. For QCD, the 
study of the infrared limit is of central importance for 
the comprehension of the mechanisms of quark and gluon 
confinement and of chiral-symmetry breaking. In what 
concerns confinement, in Landau gauge, the infrared be- 
havior of gluon andghost propagators is linked with the 
Gribov-Zwanziger [J 0] and the Kugo-Ojima W\ confine- 
ment scenarios. These confinement mechanisms predict, 
at small momenta, an enhanced ghost propagator and 
a suppression of the gluon propagator. The strong in- 
frared divergence for the ghost propagator corresponds 
to a long-range interaction in real space, which may be 
related to quark confinement. The suppression of the 
gluon propagator, which should vanish at zero momen- 
tum, implies (maximal) violation of reflection positivity 
and may be viewed as an indication of gluon conflnement. 
Moreover, the interest in the propagators goes beyond 
the confinement mechanism, as they are inputs for many 
phenomenological calculations in hadronic physics (see, 
for example, Refs. 0, §)■ 

Analytic studies of gluon and ghost propagators using 
Schwinger-Dyson equations (SDE) |S,t7|,iSi] seem to agree 
with the above scenarios. (The reader should however 
be aware that, in the literature, there are solutions of 
the SDE that do not comply with the Gribov- 

Zwanziger or the Kugo-Ojima predictions at small mo- 
menta.) Moreover, when dynamic quarks are neglected 



and assuming that ^ 1 /Nc — as suggested by analysis 
of the large Nc limit [ll| — the SDE become indepen- 
dent of the number of color Nc- Thus, they predict that 
gluon and ghost propagators be independent of Nc- 

The Landau gauge gluon propagator D{k^) has been 
investigated with lattice techniques in quenched QCD 
i.e. pure SU(3) Yang-Mills theory]_ [13, lH, H, [S [IF 
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dimensions) 
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3 and 4 space-time 
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All lattice studies 



in 4d suggest a finite nonzero infrared gluon propaga- 
tor [20|, i2^, m, 113, nil , in contradiction with the infrared 
Schwinger-Dyson solution. On the other hand, finite- 
size effects are very larg e and not yet well-controlled, 
even in the 3d case 3J|. Only in two space-time di- 
mensions [4l|, using a lattice side L up to about 40 fm, 
does one find that D{0) extrapolates to zero as L goes 
to infinity Let us note that investigation of SDE on a 
4-torus [4^ suggests that the gluon propagator indeed 
approaches the infinite-volume limit very slowly, espe- 
cially for its low-momentum components. On the other 
hand, even with an infrared-finite propagator, one clearly 
finds [13, m, HI m, SI that reflection-positivity is vio- 
lated when sufficiently large lattice volumes are consid- 
ered. Finally, in the 2d SU{2) case [4l| and in the 4d 
SU{3) case (using asymmetric lattices) [1^ [s^] it was 
found that the gluon propagator complies with the pure 
power-law behavior predicted analytically Q . 

The lattice-Landau-gauge SU{2) and SU{3) ghost 



TABLE L Lattice setup. The lattice spacing was computed 
from the string tension, assuming ^fa — 440 MeV. For SU{3), 
the lattice space was taken from 54]. The corresponding P 
values for SU (2) were computed using the asymptotic scaling 
analysis discussed in [ssj. 
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propagator Gjk"^) has been studied in |22 

M: MM, [M EE El, m, 113, El, Ei, SMI, 

in all cases an enhancement of the propagator compared 
to the tree-level behavior was observed. Concern- 
ing the comparison between lattice results and the SDE 
solution, the two propagators seem to agree only qual- 
itatively. In particular, in three and in four space-time 
dimensions, the infrared exponent obtained using lattice 
simulations is always smaller than the one predicted an- 
alytically. On the other hand, on the 2d SU{2) case [4l| . 
the ghost propagator shows an infrared behavior 
in agreement with the SDE solution Q . 

In summary, for the Landau gauge, the SDE gluon and 
ghost propagators agree, at least qualitatively, with the 
lattice propagators. However, while analytic studies us- 
ing Schwinger-Dyson equations predict the same infrared 
behavior for the SU{2) and S'C/(3) gauge groups, lattice 
simulations usually assume that the two cases are dif- 
ferent, although their qualitative infrared features may 
be the same. In this paper, we carry out a comparative 
study of lattice Landau gauge propagators for these two 
gauge groups. Our data were especially produced by con- 
sidering equivalent lattice parameters in order to allow a 
careful comparison of the two cases. We note that we 
do not assume a power-law behavior for the propagators, 
but just compare the raw data in the two cases. 



II. NUMERICAL SIMULATIONS 
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FIG. 1: Renormalized gluon propagator as a function of the 
squared magnitude of the four-momentum k, for on-axis 
momenta (fc, 0,0,0). We show results for the 5(7(3) and the 
5(7(2) cases. In the two top plots we report data for the 
three sets of lattice parameters with approximately the same 
physical lattice volume V . In the two bottom plots we report 
data for the two sets of lattice parameters with /3 = 6.0 for 
5(7(3) and (3 = 2.4469 for 5(7(2). In all cases we show only 
data for fc^ < 5 GeV'^. For larger momenta, the data using 
different lattice setups agree well. 



We consider four different sets of lattice parameters, 
with the same lattice size N'^ and the same physical lat- 
tice spacing a for the two gauge groups (see Table [l|. 
The first three cases are chosen to yield approximately 
the same physical lattice volume V = {NaY ~ (1-7 fm)^. 
This allows a comparison of discretization effects. The 
fourth case corresponds to a significantly larger physical 
volume, i.e. V ~ (3.2 fm)^, in order to study finite-size 
effects. For all four cases, 50 configurations were gen- 
erated [6lj using the Wilson action. The gluon and the 
ghost propagators 

Dtik^) - S"" ^6,. - ^) D{P) , (1) 
G"''(fc2) = -6''''G{k^) (2) 



were computed for four different types of momenta: 
(/c, 0,0,0), (fc,fc,0,0), {k,k,k,0) and {k,k,k,k). In the 
computation of D{k?) and G{k'^), an average over equiv- 
alent momenta and color components was always per- 
formed. In this work we use the field definitions and the 
choice of momenta reported in [35] for the SU{2) case 
and in 26] for SU{3). In particular, each component k 
is given (in lattice units) by A; 2sin(7r7i), where n is an 
integer. 

Here, we do not check for possible effects of the break- 
ing of rotational invariance [551 ]. In particular, we always 
compare results for the SU{2) and the 5C7(3) groups us- 
ing the same type of momenta in the two cases. We 
also do not consider possible Gribov-copy effects. In- 
deed, even though they can play an important role in the 



FIG. 2: Renormalized ghost propagator as a function of the 
squared magnitude of the four-momentum k, for on-axis 
momenta (fc, 0,0,0). The data are organized as in Fig. [1] 



FIG. 3: Ratios of 5(7(3) over SU{2) gluon propagators for 
the four lattice setups considered. 



infrared behavior of the propagators [2l|, |3l| , witli our set 
of lattice volumes and for the statistics considered here 
these effects should always be smaller than the statistical 
error. 

The propagators were computed in the minimal Lan- 
dau gauge, obtained by minimizing the functional 



s[n] 



(3) 



the statistical errors were computed with the (single- 
elimination) jackknife method in the SU{3) case and with 
the bootstrap method (using 1000 bootstrap samples) in 
the SU{2) case. We checked that these errors are in 
agreement with those obtained considering one standard 
deviation. 

In order to compare the propagators from the differ- 
ent simulations, the gluon and ghost propagators were 
renormalized accordingly to 



where U^{x) — ^l(x) Unix) ^lUx + Bn) is the gauge- n(h'^\\ — — 

J. r_^.„_j i-^i J - • J.1 1^ i 1 .. >\h^ = u^ ~ ..2' 



transformed link and is the unit vector along the /i 
direction. For SU (2) the gauge fixing was performed 
using a stochastic-overrelaxation algorithm (see [1^1 for 
details), while for SU{i) a Fourier-accelerated steepest- 
descent algorithm was used (see [26| for details). 

In what concerns the evaluation of the ghost prop- 
agator, in the SU{2) case the Faddeev-Popov matrix 
was inverted using the method described in Slf, while 
the SU{3) simulation relies on the method discussed in 
Ref. [43] (considering more than one source). In the 
calculation of the gluon and of the ghost propagators. 



G(fc2 



A' (4) 



using fi — 3 GeV as a renormalization point. The lattice 
data were interpolated (using splines) to allow the use of 
such a renormalization point in all the simulations. We 
have checked that the interpolation reproduces perfectly 
the lattice data. Let us note that, due to breaking of 
rotational invariance, the renormalization factors Z^jL^) 
depend, in general, slightly on the type of momenta. Here 
we use, for all momenta fc, the factor Z{fi^) obtained from 
the on-axis momenta (fc, 0,0,0). 
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FIG. 4: Ratios of SU{3) over SU{2) ghost propagators for 
the four lattice setups considered. 



III. RESULTS 

The renormalized SU{2) and SU{3) propagators can 
be seen for the various lattice setups in Fig. [1] (gluon) 
and in Fig. [2] (ghost) for the on-axis momenta (fc, 0,0,0). 
(Results are similar when considering the other types of 
momenta.) In all figures we report, on the horizontal 
axis, the squared magnitude (in GeV^) of the four- 
momentum k. These figures show that, for the set of 
momenta accessible in our simulations, finite- volume and 
finite-spacing effects are under control. Moreover, they 
show that the SU{2) and SU(^) propagators are essen- 
tially equal, with slight differences in the low-momenta 
region. Similar results have been recently presented at 
Lattice 2007 by Anthony G. Wilhams In Figs. E] 

(gluon) and[3] (ghost) we show the ratios of SU{3) over 



SU{2) propagators. The statistical errors were com- 
puted assuming Gaussian-error propagation. Note that 
in the case of the gluon propagator there are momenta 
for which the discrepancy from 1 for the ratio is about 
10% or larger. However, these deviations are not sys- 
tematic and are probably due to a combination of sev- 
eral effects. These may include breaking of rotational 
invariance, small statistics and finite-size effects, such as 
those related to theglobal Z(Nc) symmetry of the lattice 
action [stI. [ssI. [sqI. [6Q|] . 

IV. CONCLUSIONS 

In summary, considering a careful choice of the lattice 
parameters, we were able to carry out an unambiguous 
comparison of the lattice Landau gluon and ghost prop- 
agators for SU{2) and SU{3) gauge theories. The data 
show that the two cases have very similar finite-size and 
discretization effects. Moreover, we find very good agree- 
ment between the two Yang-Mills theories (for our values 
of momenta larger than 1 GeV), for all lattice parame- 
ters and for all types of momenta. Below 1 GeV, the 
results for the two gauge groups show some differences, 
especially for the gluon propagator. Note, however, that 
all ratios are compatible with 1 within two standard de- 
viations. 

In this sense, our results suggest that the propagators 
are the same for all SU{Nc) groups in the nonpertur- 
bative region, as predicted by Schwinger-Dyson equa- 
tions. Of course, given the lattice volumes considered, 
further studies are required before drawing final conclu- 
sions about the comparison below 1 GeV. In particular, 
it will be interesting to investigate if this agreement per- 
sists also in the deep-infrared region, where the gluon 
propagator may show a turnover and a suppression, as 
predicted in the Gribov-Zwanziger scenario. 
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